ON THE FINITE-TIME SPLASH AND SPLAT SINGULARITIES FOR 
THE 3-D FREE-SURFACE EULER EQUATIONS 
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Abstract. We prove that the 3-D free-surface incompressible Euler equations with regu- 
lar initial geometries and velocity fields have solutions which can form a finite-time "splash" 
(or "splat") singularity first introduced in [8], wherein the evolving 2-D hypersurface, the 
moving boundary of the fluid domain, self-intersects at a point (or on surface). Such 
singularities can occur when the crest of a breaking wave falls unto its trough, or in the 
study of drop impact upon liquid surfaces. Our approach is founded upon the Lagrangian 
description of the free-boundary problem, combined with a novel approximation scheme of 
a finite collection of local coordinate charts; as such we are able to analyze a rather general 
set of geometries for the evolving 2-D free-surface of the fluid. We do not assume the fluid 
is irrotational, and as such, our method can be used for a number of other fluid interface 
problems, including compressible flows, plasmas, as well as the inclusion of surface tension 
effects. 
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1. Introduction 

1.1. The Eulerian description of the free-boundary problem. For < t < T, the 

evolution of a three-dimensional incompressible fluid with a moving free-surface is modeled 
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by the incompressible Euler equations: 
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The open subset f2(t) C M 3 denotes the changing volume occupied by the fluid, F(t) := dQ(t) 
denotes the moving free-surface, V(T(t)) denotes normal velocity of T(t), and n(t) denotes 
the exterior unit normal vector to the free-surface T(t). The vector-field u = (u\, U2,us) 
denotes the Eulerian velocity field, and p denotes the pressure function. We use the notation 
D = (di, 82, 83) to denote the gradient operator. We have normalized the equations to have 
all physical constants equal to 1. 

This is a free-boundary partial differential equation to determine the velocity and pressure 
in the fluid, as well as the location and smoothness of the a priori unknown free-surface. 
In the case that the fluid is irrotational, curlu = 0, the coupled system of Euler equations 
(jl.ip can be reduced to an evolution equation for the free-surface (with potential flow in 
the interior), in which case (jl.ip simplifies to the water waves equation. We do not make 
any irrotationality assumptions in this work. 

We will prove that the 3-D Euler equations (11. 1|) admit classical solutions which evolve 
regular initial data onto a state, at finite-time T > 0, at which the free-surface self-intersects, 
and the flow map loses injectivity. The self-intersection can occur at a point, causing a 
"splash," or on a surface, creating a "splat." 

1.2. Local-in-time well-posedness. We begin with a brief history of the local-in-time 
existence theory for the free-boundary incompressible Euler equations. For the irrotational 
case of the water waves problem, and for 2-D fluids (and hence 1-D interfaces), the earliest 
local existence results were obtained by Nalimov [21] , Yosihara [32] , and Craig [10] for initial 
data near equilibrium. Beale, Hou, & Lowengrub [5] proved that the linearization of the 
2-D water wave problem is well-posed if the Taylor sign condition 



is satisfied by the initial data (see [23] and (25]). Wu [28] established local well-posedness 
for the 2-D water waves problem and showed that, due to irrotationality, the Taylor sign 
condition is satisfied. Later Ambrose & Masmoudi [3], proved local well-posedness of the 
2-D water waves problem as the limit of zero surface tension. For 3-D fluids (and 2-D 
interfaces) , Wu }29] used Clifford analysis to prove local existence of the water waves problem 
with infinite depth, again showing that the Taylor sign condition is always satisfied in the 
irrotational case by virtue of the maximum principle holding for the potential flow. Lannes 
|19j provided a proof for the finite depth case with varying bottom. In all of these works, 
the free-surface was considered a perturbation of a flat interface. 



dp 
dn 
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The first local well-posedness result for the 3-D incompressible Euler equations without 
the irrotationality assumption was obtained by Lindblad [20] in the case that the domain 
is diffeomorphic to the unit ball using a Nash-Moser iteration. In Coutand and Shkoller 
[13] . we obtained the local well-posedness result for arbitrary geometries that have H - 
class boundaries and without derivative loss (this framework, employing local coordinate 
charts in the Lagrangian configuration, is ideally suited for the splash and splat singularity 
problems that we study herein). Shatah and Zeng [25] established a priori estimates for 
this problem using an infinite-dimensional geometric formulation, and Zhang and Zhang 
proved well-poseness by extending the complex-analytic method of Wu [29] to allow for 
vorticity. Again, in the latter case the geometry was constrained to be a perturbation of 
the half-space. 

1.3. Long-time existence. It is of great interest to understand if solutions to the Euler 
equations can be extended for all time when the data is sufficiently smooth and small, or if 
a finite-time singularity can be predicted for other types of initial conditions. 

Because of irrotationality, the water waves problem does not suffer from vorticity con- 
centration; therefore, singularity formation involves only the loss of regularity of the inter- 
face. In the case that the irrotational fluid is infinite in the horizontal directions, certain 
dispersive-type properties can be made use of. For sufficiently smooth and small data, 
Alvarez- Samaniego and Lannes [2] proved existence of solutions to the water waves prob- 
lem on large time-intervals (larger than predicted by energy estimates), and provided a 
rigorous justification for a variety of asymptotic regimes. By constructing a transforma- 
tion to remove the quadratic nonlinearity, combined with decay estimates for the linearized 
problem (on the infinite half-space domain), Wu |30] established an almost global existence 
result (existence on time intervals which are exponential in the size of the data) for the 2-D 
water waves problem with sufficiently small data. Wu |31] then proved global existence in 
3-D for small data. Using the method of spacetime resonances, Germain, Masmoudi, and 
Shatah [T7] also established global existence for the 3-D irrotational problem for sufficiently 
small data. 

1.4. Splasing of liquids and the finite-time splash singularity. The study of splash- 
ing, and in particular, of drop impact on liquid surfaces has a long history that goes back to 
the end of the last century when Worthington |27] studied the process by means of single- 
flash photography. Numerical studies show both fascinating and unexpected fluid behavior 
during the splashing process (see, for example, Oguz & Prosperctti [22J), with agreement 
from matched asymptotic analysis by Howison, Ockendon, Oliver, Purvis and Smith [18]. 

The problem of rigorously establishing a finite-time singularity for the fluid interface has 
recently been explored for the 2-D water waves equations by Castro, Cordoba, Fefferman, 
Gancedo, and Gomez-Serrano in [U [9], where it was shown that a smooth initial curve 
exhibits a finite-time singularity via self-intersection at a point; they refer to this type of 
singularity as a "splash" singularity, and we will continue to use this terminology. (We will 
give a precise definition of the splash domain in our 3-D framework in Section 13.1.21 and we 
define the splat domain in Section [9j) 

Their work follows earlier results by Castro, Cordoba, Fefferman, Gancedo, and Lopez- 
Fernandez [6] and Castro, Cordoba, Fefferman, Gancedo, and Lopez-Fernandez [7] for both 
the Muskat and water waves equations wherein the authors proved that an initial curve 
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Figure 1. The splash singularity wherein the top of the crest touches the 
trough at a point xq in finite time T. 



which is graph, that satisfies the Taylor sign condition, reaches a regime in finite time in 
which it is no longer a graph and can become unstable due to a reversal of the sign in the 
Taylor condition. 

Herein, we develop a new framework for analyzing the finite-time splash and splat sin- 
gularity for 3-D incompressible fluid flows with vorticity. Our motivation is to produce a 
general methodology which can also be applied to compressible fluids, as well as to ionized 
fluids, governed by the equations of magnetohydrodynamics. Our method is founded upon 
the transformation of (II. ip into Lagrangian variables. We are thus not restricted to poten- 
tial flows, nor to any special geometries. Furthermore, our method of analysis does not, in 
any significant way, distinguish between flow in different dimensions. While we present our 
results for the case of 3-D fluid flow, they are equally valid in the 2-D case. 



1.5. Main result. The main result of this paper states that there exist initial domains Qq 
of Sobolev class H 4 5 together with initial velocity vectors no £ H 4 (Uq) which satisfy the 
Taylor sign condition f 1 1 . 2 1> . such that after a a finite-time T > the solution of the Euler 
equations reaches a "splash" (or "splat") singularity. At such a time T, particles which 
were separated at time t = collide at a point xq (or on a surface To), the flow map r/(T) 
loses injectivity, and d[£l c ] forms a cusp. In short, T is the time at which the crest of a 3-D 
wave turns-over and touches the trough. This statement is made precise in Theorems 15.11 
and EH 

Note that the use of i? 4 5 -regularity for the domain Qq and -ff 4 (f2o)- re g u l ar ity for velocity 
field uq is due to the functional framework that we employ for the a priori estimates in 
Theorem IA.1I For 3-D incompressible fluid flow, we find that this is the most natural 
functional setting; of course, we could also employ any .fP-framework for s > 4.5 or a 
Holder space framework as well. 



1.6. The Lagrangian description. We transform the system (jl.l|) into Lagrangian vari- 
ables. We let rj(x,t) denote the "position" of the fluid particle x at time t. Thus, 



dtf] = u o r] for t > and rj(x, 0) = x 
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where o denotes composition so that [u o r]](x,t) := u(n(x,t),t) . We set 

v = u o rj (Lagrangian velocity), 
q = p o r] (Lagrangian pressure) , 
A = [Dt]]" 1 (inverse of the deformation tensor), 
J = det[Dr]] (Jacobian determinant of the deformation tensor), 
a = J A (cofactor of the deformation tensor). 

Whenever divu = 0, it follows that detDr] = 1, and hence the cofactor matrix of Drj is 
equal to [Drj\~ l , i.e., a = A. Using Einstein's summation convention, and using the notation 
F,k to denote J^, the ki/i-partial derivative of F for k = 1,2,3, the Lagrangian version of 
equations (jl.ip is given on the fixed reference domain fi by 



v (t) = e+[ v inOx [0,T], (1.3a) 
Jo 

v t + A T Dq = infix (0,7*], (1.3b) 

div v v = in x [0,7] , (1.3c) 

q = onrx[0,T], (1.3d) 

(rj,v) = (e,u ) in fi x {t = 0} , (1.3e) 



where e(x) = x denotes the identity map on fi, and where the ith-component of A T Dq is 
A\ (A T denotes the transpose of the matrix A.) By definition of the Lagrangian flow 
rj(t), the free-surface is given by 

r(t) = f|(t)(r) . 

We will also use the notation rj(t, T) = T(t), and rj(t, fi) = fi(t). The Lagrangian divergence 
is defined by div v v = A\v l ,j. Solutions to (II. 3p which are sufficiently smooth to ensure 
that r](t) are diffeomorphisms, give solutions to (jl.lj) via the change of variables indicated 
above. 

1.7. The splash singularity for other hyperbolic PDEs. Our methodology can be 
applied to a host of other time-reversible PDEs that have a local well-posedness theorem. 

(1) Surface tension. Our main result also holds if surface tension is added to the Euler 
equations. In this case equation (jl.3d|) is replaced with 

qn = -aA g (rj) , 

where a > denotes the surface tension parameter, n is the outward unit-normal 
to T(t), A g denotes the surface Laplacian with respect to the induced metric g 
where g a p = rj, a -r],p. This is the Lagrangian version of the so-called Laplace- Young 
boundary condition for pressure: p = o~H, where H is the mean curvature of the 
free-surface T(t). We have established well-posedness for this case in [13]. The only 
modifications required for the case of positive surface tension is to consider initial 
domains fio of Sobolev class H 6 with initial velocity fields uq G i7 4,5 (fio)- Our main 
theorem then provides for a finite-time splash singularity for the case that a > 0. 
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(2) Physical vacuum boundary of a compressible gas. We can also consider the evolution 
of the free-surface compressible Euler equations which model the expansion of a gas 
into vacuum. We established the well-posedness of this system of degenerate and 
characteristic multi-D conservation laws in [15] . In this setting, we our methodology 
shows that there exist initial domains Qq of class H 4 , initial velocity fields uq £ 
ff 3 ' 5 (f2o)> an d initial density functions po £ H 4 (Qq) such after time T > 0, a splash 
singularity if formed by the evolving vacuum interface. 

(3) Other physical models. In fact, we can establish existence of a finite-time splash sin- 
gularity for a wide class of hyperbolic systems of PDE which evolve a free-boundary 
in a sufficiently smooth functional framework, and which are locally well-posedness. 
Examples of equations (not mentioned above) include nonlinear elasticity and mag- 
netohydodynamics . 



used to denote the three-dimensional gradient vector D = i - . I 



2. Notation, local coordinates, and some preliminary results 

2.1. Notation for the gradient vector. Throughout the paper the symbol D will be 

d d d \ 
dxi ' 8x2 ' 8x3 J ' 

2.2. Notation for partial differentiation and Einstein's summation convention. 

The kih partial derivative of F will be denoted by F,^ = J^-. Repeated Latin indices i, j, k, 
etc., are summed from 1 to 3, and repeated Greek indices a,/3,j, etc., are summed from 1 
to 2. For example, F, u = £* = i agfc and F\ a I^G\, = ELi ELi E?=i f£^fg- 

2.3. The divergence and curl operators. For a vector field u on O, we set 

divu = u 1 ,! +u 2 , 2 +u 3 , 3 , 
curln= (u 3 , 2 -u 2 , 3 ,u 1 , 3 -u 3 , 1 ,v?,i-u 1 ,2) ■ 

[ 1, even permutation of {1, 2, 3}, 

With the permutation symbol eijk given by e ijk = I -1, odd permutation of {1,2,3}, the ith- 

0, otherwise, 

component of curln is given by 

(curltt)j = e ijk u k ,j . 



2.4. The Lagrangian divergence and curl operators. We will write div^ v = divuorj 
and curl^ v = curl uo r/. From the chain rule, 

drvv, v = A s r v r , s and (curl^ v)i = Eij h A S jV k , s . 

2.5. Local coordinates near T. In Appendix A, we establish the a priori estimates for 
solutions of the 3-D free-surface Euler equations (following our local well-posedness theory 
in |131ll4j). Such solutions evolve a moving two-dimensional surface which is of Sobolev class 
H 4 . This boundary regularity implies a three-dimensional domain of class i? 4,5 , constructed 
via a collection of i/ 4,5 -class local coordinates. 
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Let Hcl 3 denote an open subset of class and let {Ui}fL 1 denote an open covering 
of r = d£l, such that for each I £ {1, 2, . . . , K}, with 

B = B(0, 1), denoting the open ball of radius 1 centered at the origin and, 
B + = Bn{x 3 > 0}, 
B° = Bn{x 3 = 0}, 

there exist iif 45 -class charts Q[ which satisfy 

9f. B — > XJ\ is an if 4,5 diffeomorphism, (2.1a) 

e l (B + ) = u l nn, 9 l (B°) = u l nv. (2.1b) 

Next, for L > K, we let {U{\f =K+1 denote a family of open sets contained in Q such that 



( fK/f ), 

Uk\ "X ) ''K-l ■ • ■ ] 

V -Vt^rnt 



Figure 2. Indexing convention for the open cover {U{\f =l of 0. 

1 is an open cover of Vt, and we such that there exist diffeomorphisms Oi : B —> Ui. 

2.6. Tangential (or horizontal) derivatives. On each boundary chart Ui H f2, for 1 < 

I < K, we let d denote the tangential derivative whose ath-component given by 

For functions defined directly on B + , d is simply the horizontal derivative d = (d Xl ,d X2 ). 

2.7. Sobolev spaces. For integers k > and a domain E of M 3 , we define the Sobolev 
space H k (E) (H k (E;R 3 )) to be the completion of C°°(E) (C°°(E; ]R 3 )) in the norm 

\w\\Ie= E / i^ a ^)i 2 > 

\a\<k JE 

for a multi-index a G Z 3 , , with the convention that |a| = ai + a2 + (23. When there is no 
possibility for confusion, we write || • \\k for || • \\kE- For real numbers s > 0, the Sobolev 
spaces H S (E) and the norms || • \\ Sj e are defined by interpolation. We will write H S (E) 
instead of H S (E;M. 3 ) for vector-valued functions. 
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2.8. Sobolev spaces on a surface T. For functions u £ H k 



(r), k > 0, we set 



<v = E / \ Ba < x ) 



\a\<k 



for a multi-index a £ Z+. For real s > 0, the Hilbert space H S (T) and the boundary norm 
| • \ s is defined by interpolation. The negative-order Sobolev spaces H~ S (T) are defined via 
duality. That is, for real s > 0, H~ S (T) = H S (T)'. 

2.9. The norm of a standard domain fL 

Definition 2.1. ^4 domain Q, is of class H 45 if for each I = 1, ...,L, each diffeomorphism 
0i is of class H 45 . The H 4,5 -norm of Q, is defined by 



In particular if e : £1 — > ft is the identity map, then He^.s^ is given by 

We can, of course, replace H 4 5 with any H s , s > 2.5 to define domains of class H s . 

2.10. Local well-posedness for the free-surface Euler problem. 

Theorem 2.1 (Coutand and Shkoller [13]). With E(t) given by jlA.8\) . suppose that E(0) < 
Mq and that the initial pressure function satisfies the Taylor sign condition. Then there 
exists a solution to ( fi. 1\) on [0,T] where T > depends E(0), and sup t6 [ r ] E(t) < 2Mq. 

3. The splash domain ft s and its approximation by standard domains $7 e 
3.1. The splash domain. 

3.1.1. The meta- definition. A splash domain f2 s is an open and bounded subset of M. n which 
is locally on one side of its boundary, except at a point xq £ d£l s , where the domain is locally 
on each side of the tangent plane at xq. The domain Q s satisfies the cone property and 
can be approximated (in sense to be made precise below) by domains which have a smooth 
boundary. 

We observe that the Sobolev spaces H r (Q s ) are defined for the splash domain £l s in the 
same way as for a domain which is locally on one side of its boundary; moreoever, as the 
bounded splash domain Q s satisfies the cone property, interpolation theorems and most 
of the imporant Sobolev embedding results hold (see, for examples, Chapters 4 and 5 of 
Adams pQ). 

The main difference between bounded splash domains with the cone property and do- 
mains that have the uniform i? r -regularity property is with regards to trace theorems: For 
the splash domain Q, s , a function / in H 4 ' 5 (Q S ) has a trace in H 4 (T') for any smooth subset 
r' of d£l s whose closure does not contain xq. At xo there is not a well-defined (global) 
trace for /, in the sense of coming from both sides of the tangent plane at xq, although it is 
indeed possible to define local traces for / at xq with respect to each of the local coordinate 
charts containing xq. 




(2.2) 
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3.1.2. The definition of the splash domain. 

(1) We suppose that xq G V := d£l s is the unique boundary self-intersection point, 
i.e., Q s is locally on each side of the tangent plane to d£l s = T s at xq. For all 
other boundary points, the domain is locally on one side of its boundary. Without 
loss of generality, we suppose that the tangent plane at xq is the horizontal plane 
^3 - (2:0)3 = 0. 

(2) We let Uq denote an open neighborhood of xq in R 3 , and then choose an additional 
L open sets {Ui}f =l such that the collection {Ui}f =Q is an open cover of T s , and 

is an open cover of Q s and such that there exists a sufficiently small open 
subset uj C Uq containing xq with the property that 



u+ = u nn s r\ {x 3 > (x ) 3 } and u =u nn s r\ {x 3 < {x ) 3 } . 

Additionally, we assume that Uq n O s n {2:3 = (2:0)3} = {2:0}, which implies in 
particular that Uq and L^q - are connected. 



Figure 3. Splash domain Q s , and the collection of open set 
{U , Ui, U 2 , U K } covering T. 

(3) For each / G {1, ...,K}, there exists an i7 4 ' 5 -class diffeomorphism 61 satisfying 



mnUi = Hl for all l = l,...,L. 



We set 




6r.B:= B(0, 1) -> Ui 
6i{B+) and n T s = 0,(5°) , 



where 



B+ 



{(xi,x 2 ,x 3 ) G B : x 3 > 0} , 
{(xi,x 2 ,x 3 ) G 5 : x 3 = 0} . 



Z? 1 
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(4) For L > K, let {U{\f =K+1 denote a family of open sets contained in fi s such that 
{U{\i = q is an open cover of Q s , and for I £ {K + 1, ...,L}, 9\ : B — > U\ is an i? 45 
diffeormorphism. 

(5) To the open set f7o we associate two H 5 -class diffeomorphisms 9 + and 0_ of B 
onto [/o with the following properties: 

6 + (B + ) = U+, 9^(B + ) = U^, 

9 + (B°) = u+nr s , d^(B°) = u^nr s , 

such that 

{x o } = e + (B°)n0-(B ), 

and 

e + (o) = e^(o) = x . 

We further assume that 

6±(B+ n 5(0, 1/2)) n 6i(B+) = for Z = 1, 

and 

0±(.B + n B{0, 1/2)) n 0i(5) = for I = K + 1, L . 

Definition 3.1 (Splash domain J7 S ). We say that £l s is a splash domain, if it is defined 
by a collection of open covers {Ui}(l and associated maps 9i, 02, Ol,} satisfying the 
properties (l)-(5) above. Because each of the maps is an H ° diffeomorphism, we say that 
the splash domain O s defines a self-intersecting generalized H 45 -domain. 

3.2. A sequence of standard domains approximating the splash domain. We ap- 
proximate the two distinguished charts 0- and 9+ by charts 01 and 01 in such a way as to 
ensure that 

9t{B°)r\9 e + (B°) = V e > 0, 

and which satisfy 

9 e _ 9- and 9% ->■ 9+ as e ->• . 
We choose r > sufficiently small so that 

0_(B + (O,2t)) c u and 0+(B + (0,2v)) Cu, 

and then we let ip £ T>(B(0,v)) denote a smooth bump-function satisfying < tp < 1 and 
ip(0) = 1. For e > taken small enough, we define 

9l(x) = 0_(x) -e V(x)e 3 , 
0^(x) = 0+(x) + e ^(x)e 3 , 

where e 3 = (0,0, 1) denotes the vertical basis vector of the standard basis e; of M 3 . By 
choosing ip £ T>(B(0,z)), we ensure that the modification of the domain is localized to a 
small neighborhood of xq and away from the boundary of Uq and the image of the other 
maps 9\. Then, for e > sufficiently small, 

9t(B+) n 9 e + (B+) = 0. 




e! 

Figure 4. The black dot denotes the point xq where the boundary self- 
intersects (middle). For e > 0, the approximate domain £l e does not intersect 
itself (right). 



Since the maps 9± are a modification of the maps 0± in a very small neighborhood of G B, 
we have that for e > sufficiently small, 

9 e ± {B + n B(0, 1/2)) n 9 t (B + ) = for I = 1, K , 

and 

0±(5 + n B(0, 1/2)) n 9i(B) = for J = K + 1, L . 

For Z G {1, we set 0f = 9 h Then (91 : B+ -> C/ , 0+ : # + -> #0, and 0f : 5+ -> £/j, 

/ € {1, ...,-ftT}, Of : B Ui, I G {-ftT + 1, ...,L}, is a collection of H 4 5 coordinate charts as 
given in Section 12.51 and so we have the following 

Lemma 3.1 (The approximate domains £l e ). For each e > sufficiently small, the set fi e , 
defined by the local charts 6*1 : B + -> U , 9% : B + ->■ f/ , a™d 0f : B + ^ U h I £ {1, K}, 
Of : B —> Ui, I G {-ftT + 1, L} is a domain of class H 4 " 5 , which is locally on one side of its 
H boundary. 

By choosing < ro < \ such that 1 > tp > \ in B(0, ro), we see that 

| (01 (a?) - 0+(y)) • e 3 | > e for any x, y G 5+ n 5(0, r ) . 

With ro chosen, due to the fact that by assumption (2) the images of 0_ and + only intersect 
the plane {X3 = (^0)3} at the point xq, there exists 5(ro) > such that (9 e (x) — xq) • e% < 
—5(ro) and (0 e (x) — xq) • e 3 > 6(ro) for all x G -B + with |x| > ro. This, in turn, implies that 
if x G B + with I a; I > ro and y G B + , we then have that 

|(01(x)-0;(y))-e 3 |> ( 5(ro)-26>e if e < ^ . 

We have therefore established the following fundamental inequality: for < e < ^g * 1 , 

V(x,y) e B+ x B+ , \(9t{x)-9Uy))-e 3 \>e. (3.1) 
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We henceforth assume that < e < g , 

In summary, we have approximated the self-intersecting splash domain £l s with a sequence 
of i? 4 - 5 -class domains f2 e , < e < (such that dQ e does not self- intersect). As such, each 
one of these domains f2 e , e > 0, will thus be amenable to our local-in-time well-posedness 
theory for free-boundary incompressible Euler equations with Taylor sign condition satisfied. 

We also note that f2 e and £l s are the same domain, except on the two patches 9 e _{B + n 
B(0, |)) and 9+(B + nB(0, |)). In particular, as 8± differ from 0± on a set properly contained 
in to C Uq, we may use the same covering {Ui}fL Q for Q e as for Q s . 

Lemma 3.2. For < e < ^ ' , i/ie B 4 ^ -norm of Q e is bounded independently of e. 

Proof. The assertion follows from ||0±||4.5.b+ < ||0±||4.5,b+ + jj ||V'I|4.5,.B+ • ^ 

3.3. A uniform cut-off function on the unit-ball B. Let B\^ a = B(0, 1 — a) for < 
a < 1. For a > taken sufficiently small, we have that 9-(Bf_ a ) C Co and 9+(Bf__ a ) C £/o 
and for each / = 1, K, 9[{Bf_ a ) C Ui, and for each I = K + 1, L, 8i{B\_ a ) C C/j, and 
the open sets 8„(B+_ a ), 9 + (B+_ a ), 9 t (B+_ a ) (1 < Z < K), tf^i-a) (JiC + 1 < Z < L), are 
also an open cover of f2 s . Since the diffeomorphisms 0|. are modifications for 9± in a very 
small neighborhood of the origin, it is clear that independently of e > 0, the sets 9L(Bf_ a ), 
9 +( B i~a)> °i( B t~a) (1 < ' < K ), 0l{Bi-a) (K + 1 < I < L) are also an open cover of each 

Definition 3.2 (Uniform cut-off function (). Let ( G £>(-B(0, 1)) such that < ( < 1 and 
((x) = 1 for \x\ < 1 — a and ( = /or |x| > 1 — 
We set c = 1 — so i/iai 

0<C€2?(J3(0,0)<1- (3-2) 

4. Construction of the splash velocity field u s at the time of the splash 

singularity 

We can now define the so-called splash velocity u s associated with the generalized H 4 5 - 
class splash domain Q s , as well as a sequence of approximations u e s set on our i? 4,5 -class 
approximations O e of the splash domain £l s . 

4.1. The splash velocity u s . 

Definition 4.1 (Splash velocity u s ). A velocity field u s on an H 4 ' 5 -class splash domain Q s 
is called a splash velocity if it satisfies the following properties: 

(1) (u s o 9± G B 4 - 5 (B + ), (u s o 0, G H 4 - 5 (B+) for each 1 < I < K and u s G B 4 - 5 (uj) for 
each U C f2 s ; 

so that under the motion of the fluid, the sets Uq and U ~ relatively move towards 
each other, we require that 

ujj o 0_ > CL , -u 3 s o9 + >C + in B + and C_ + C+ > , (4.1) 

where C_ and C+ are constants. 
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Definition 4.2 (Splash pressure p s ). A pressure function p s on an H 4 5 -class splash domain 
fl s is called a splash pressure associated to the splash velocity u s if it satisfies the following 
properties: 

(1) p s is the unique solution of 

dul du 3 s „ , , n . 

-Ap s = - s m n 8 , 4.2a 

OXj OXi 

Qp s o e± = on B° , (4.2b) 

Cp,o0, = O onB° for l = l,...,K; (4.2c) 

(2) the splash pressure p s £ H 4 - 5 (Q S ) and satisfies the local version of the Taylor sign 
condition: 

— (C Ps o8±)>0 and — ((ps o9 l )>0 on B° for I = 1, K . (4.3) 
0x3 0x3 

Note that the outward unit normal to dB + n B° points in the direction of —e^. 

Remark 1. As x = 0-(0) = 0+(0), and as p o 6L(0) = po 6 + (0) = 0, the conditions 
(4-2b\ ) and {4-% c \j are equivalent to having the usual vanishing trace p = on T s - As such, 



pe^(n s )nff i5 (S] s ). 

For property (1) in Definition 14.21 we note that p s is the unique Hq(Q, s ) weak solution 
of (|4.2p guaranteed by the Lax-Milgram theorem in £l s - The usual methods of elliptic 
regularity theory show that CPs 0± aud each CPs £ H^(B + ) for I = 1, .., L, and thus 
that p s € H 4 ' 5 (£l s )- (Notice that it is the regularity of our charts 6± and 9\ which limits 
the regularity of the splash pressure p s .) 

As we have defined in property (2) of Definition 14.21 a t the point of self-intersection xq, 
the gradient Dp s has to be defined from each side of the tangent plane at xq] namely, we 
can define Dp s o 9- and Dp s o 9 + on B°, and these two vectors are not equal at the origin 
which is the pre- image of xq under both 9- and 9+. 

It is always possible to choose a splash velocity u s so that (|4.3p holds. For example, if 
we choose u s to satisfy curlu s = 0, then (|4.3p holds according to the maximum principle 
[281 129j . On the other hand, it is not necessary to choose an irrotational splash velocity, 
and we will not impose such a constraint. Essentially, as long as the velocity field induces 
a positive pressure function, then (|4.3p is satisfied. 

4.2. A sequence of approximations u e s to the splash velocity. For e > 0, we proceed 
to construct a sequence of approximations u e s : Q e — > M 3 to the velocity field u s : O s — > M 3 
in the following way: 

u e s o 81 = u s o 9i , in B + , for I = 1, K ; (4.4a) 
u e s o 81 = u s o 0i , in B , for I = K + 1, ...,L ; (4.4b) 
u\ o 9t = u s o 9- , and , u\ o 9\ = u s o 9 + , in B + . (4.4c) 
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We then have the existence of constants A > 0, B > such that 
Klk5,n« < A (1104 o 0L|| 4|B+ + ||C< o #;|| 4 ,b+ 

K L \ 

+ EH^°^ll4.5,B + + E HC<°^ll4.5,B ] <S||u a ||4.5,n.- (4-5) 

z=i i=i<r+i / 

We next define the approximate pressure function p e s in Cl e as the ^(ST) weak solution of 

= on dn e . (4.6b) 

Again, standard elliptic regularity theory then shows that p e s E (ST). Furthermore, 
since 9± — > 9± and Of — > 9i in H 4 ^(B + ), we infer from the definition of u e s in (|4.4p that 
o 9| -> (p o fl ± and o 0f — > (p o 6i in H A ^(B + ). We may thus conclude from the 
pressure condition (|4.3|) that we also have, uniformly in e > small enough, that 

■^—(CpI ° 0±) > and — — ((p e s o Of) > , on 5° for each 1<1<K. (4.7) 

4.3. Solving the Euler equations backwards-in-time from the final states Q e and 

u e s . Because the Euler equations are time-reversible, we can solve the following system of 
free-boundary Euler equations backward-in-time: 

rf(t) =e+ f v e in fl € x [-T £ ,0] , (4.8a) 
Jo 

v e t + [A e ] T Dq £ = in O e x [-T £ , 0) , (4.8b) 

drvy v £ = in O e x [-T e , 0] , (4.8c) 

g e = on r £ x[-T £ ,0], (4.8d) 

(r] e ,v e ) = (e,u e s ) infi £ x{t = 0}, (4.8e) 

where ^4 e (x,i) = \Drf(x, t)]~ l . Thanks to Lemma 19.11 (14. 5p . and (14. 7p . we may apply 
our local well-posedness Theorem 12.11 for (14. 8|) backward-in-time. This then gives us the 
existence of T e > 0, such that there exists a Lagrangian velocity field 

v e E L°°(-T £ ,0;tf 4 (O £ )), (4.9) 

and a Lagrangian flow map 

E L°°(-T e , 0; tf 4 - 5 (JT))) (4.10) 

which solve the free-boundary Euler equations (|4.8p with /ma/ data ti £ and /inaZ domain Sl e . 
Denoting the corresponding Eulerian velocity field by 

u e = v e orf- 1 , (4.11) 

it follows that u £ is in L°°((— T e , 0); i? 4 (f2 e (t))), where f2 e (£) denotes the image of £l £ under 
the flow map rj e (t). 

In the remainder of the paper we will prove that the time of existence T € > (for our 
sequence of backwards-in-time Euler equations) is, in fact, independent of e; that is, T e 
is equal to a time T > 0, and that ||u £ (i) ||//4(n e (t)) an d ||^ e (^)||_ff 4 - 5 (r2 e ) are bounded on 
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[— T, 0] independently of e. This will then provide us with the existence of a solution which 
culminates in the splash singularity Q s at t = 0, from the initial data 

uo = lim u e (—T) , 

£->0 

Q = lim Q e (-T). 

e->0 

In particular, when solving the Euler equations forward- in-time from the initial states Qq 
and uo, the smooth H 4 ' 5 domain ^o is dynamically mapped onto the H 4,5 -class splash do- 
main £l s after a time T, and the boundary "splashes onto itself" creating the self-intersecting 
splash singularity at the point xq. 



5. The main results 

Theorem 5.1 (Finite-time splash singularity). There exist initial domains of class if 4,5 
and initial velocity fields uq E H 4 (Qq), which satisfy the Taylor sign condition \1.2) . such 
that after a finite time T > 0, the solution to the Euler equation rj(t) (with such data) maps 
rio onto the splash domain S; satisfying Definition \3.1l with final velocity u s . This final 
velocity u s satisfies the local Taylor sign condition on the splash domain Q, s in the sense of 
$i4-3\ )- The splash velocity u s has a specified relative velocity on the boundary of the splash 
domain given by {4-lty - 

The proof of Theorem [5] is given in Sections [6HHJ In Sections [9H10I we define the splat 
domain Q s and associated splat velocity u s and establish the following 

Theorem 5.2 (Finite-time splat singularity). There exist initial domains Oo of class H 4 5 
and initial velocity fields uq E H 4 (Qq), which satisfy the Taylor sign condition \1.2) . such 
that after a finite time T > 0, the solution to the Euler equation rj(t) (with such data) maps 
0,q onto the splat domain fi s , satisfying Definition \9.1\ with final velocity u s . This final 
splat velocity u s satisfies the local Taylor sign condition on the splat domain fi s in the sense 
°f i4-3\ )- The splat velocity u s has a specified relative velocity on the boundary of the splat 
domain as stated in Definition \10.1\ 



6. Euler equations set on a finite number of local charts 

For each e > 0, the functions v e , r] e , and u e , given by (|4.9p - (|4.11|) . are solutions to the 
Euler equations (|4.8p on the time interval [— T e ,0]. 

For the purpose of obtaining estimates for this sequence of solutions which do not depend 
on e > 0, we pull-back the Euler equations (|4.8p set on ST by our charts Q\. and Of, 
I = 1, L; in this way we can analyze the equations on the half-ball B + . 

It is convenient to extend the index I to include both I = and I = —1; in particular, we 
set 

01i = 0- and 0_i = 6L , 
#o = 0+ 0o = 0+ • 
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Furthermore, since for I = —1,0, 1,2, ...,K, the domain of 9\ is the half-ball B + , and for 
I = K + 1, L, the domain of 9\ is the unit-ball B, it is convenient to write 

6^ : B ^ U , 6+ : B ^ , 6 r . B ^ Ui HQ, for 1 = 1, K , 

Or.B-tUi for I = K + 1,...,L, 

so that B denotes B + for I = —1, 0, 1, 2, if and B denotes B for Z = K + 1, L. 
The Euler equations, set on B, then take the following form: 

v e = e+[ v e in Q e x [-T e ,0) , (6.1a) 

d t u e o 0f + [b\] T D{q e o 0f ) = in B x [-T e , 0) , (6.1b) 

div^oef v e o Of = in B x [-T e , 0) , (6.1c) 

q e o6f = on .Bq x [-T e , 0) , (6.1d) 

(?f o 91 v e o 0f) = (0f,u e a o 0f) onBx{t = 0}, (6.1e) 

where denotes the transpose of the matrix b^, and where for any I = —1,0, 1,2, ...L, 
M(x,t) = [D(r] e (9f(x),t)]~ . For I = K + 1,...,L, the boundary condition (|6.1d|) is not 
imposed. 

The system (|6.1|) will allow us to analyze the behavior of r) e , v e , and q € in an e- independent 
fashion. Fundamental to this analysis is the following 

Lemma 6.1 (Equivalence-of-norms lemma). With the smooth cut-off function £ given in 
Definition \3.SX there exist constants C\ > and C2 > such that for any e > and 
f £ H s (tl) with < s < 4.5, 

L L 

Ci £ IK/ flffo < Il/ll2,n« lie/ o 0f ||? iB . (6.2) 

i=-i i=-i 

Proof. Since by construction ||#f|| 4.5,3 < Cj, the first inequality is obvious. For the second 
inequality, we simply notice that with E = {x G B\ ((x) = 1}, Q e = {jfL_ 1 9f(E), so that 

L L 

11/L.n- < C J] ||/IU W < C J] ||C((^)- 1 )/II^(E) 
z=-i l=-l 
L L 

<cY. IIC/^DIU^IK^)- 1 !!^^) < c \\cm)h,B 

l=~l l=-X 
where we used the fact that detD0[ > q > for the last inequality. □ 

7. Time of existence — T of solutions to ( j4.8[ ) is independent of e 

Recall that for e > 0, the functions v e , rj £ , and u e , given by (I4.9p - (l4.1ip . are solutions 
to the Euler equations (|4.8j) on the time interval [— T e ,0]. We now prove that the time of 
existence — T € is, in fact, independent of e. 
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We begin by using the fundamental theorem of calculus to express the difference between 
the flow of two particles x and y as 

rj e (x,t) - r] e (y,t) = x - y + / [v e (x, s) - v e (y, s))ds . 



Jo 

Next, for any x and y in Q £ for which we do not have at the same time x S 9 e _(B + ) and 
y £ 8^(B + ), we see that independently of e > small enough, 

\r ] e (x,t)-r ] t (y,t)-{x-y)\<C 1 \t\ sup E e {t) \x - y\ , (7.1) 

[-T',0] 

where we have used the Sobolev embeddding theorem and where 

E e (t) = ||T/ 6 (t)||i 5 ,n« + \\v e (t)\\i,n* + ||curlu e (t)||§. 5i n« + ||«f(*)||!.5,n« ■ 

The inequality (|7.ip cannot be independent of e > if both x G 6'L(B + ) and y 6 9 e + (B + ), 
for in this case, according to (|3.ip . \x— y\ = 0(e) as e — > 0, whereas |f e (x,t) — = 0(1) 

as e — )• 0, and this, in turn, yields a global Lipschitz constant for v € of O(^) as e — >• 0. 

When x S e _(B + ) and y E #^(1? + ), there exist constants C_, C + , and a polynomial 
function Pi which are each independent of e, such that 

\ V € {x,t)- V € {y,t)\ > \ { V € {x,t)- V *{y,t))-e 3 \ 

> \ (x - y) ■ e 3 + t[u e s (x) - u e s (y)] ■ e 3 | 

-|e 3 - / v € (x,t')-u e s (x)dt'\-\e 3 - [ v e {y,t') - u e s {y)dt'\ 



o 



Xi x 2 

+ 2; 



> e + (C_ + C+)|t| - t^Pi( sup P e ) , (7.2) 

[-T E ,0] 

where the triangle inequality has been employed together with (|4.ip and (|3.ip . In order to 
obtain the lower bound on the terms X\ and I2, we again use the fundamental theorem of 
calculus, and write 

v e (x,t') = u\(x) + / v e t (x,T)dT = u e s (x) - [A e (x,T)] T Dq E (x,r)dT ; 
Jo Jo 

using (|4.8bp . it follows that 

||^(-,t')-<(-)llL-(^) < / WA'hMDfh^dr <Ct'P 1 { sup P e ) . 

JO h^ E ,o] 

We proceed to show how the two inequalities (|7.ip and (|7.2|) (together with the fact that 
C_ > and C + > 0) are used to prove that the time — T is independent of e, the flow map 
rf is injective on [— T, 0], and the a priori estimates for solutions of (|4.8p are independent 
of e on [-T,0]. 

We first record our basic polynomial-type a priori estimate, given in Theorem lA.il in the 
appendix (see also P3][I1]); we find that on [— T e ,0], 

sup E e (t) < M e + tP 2 ( sup E e (t)) , (7.3) 
te[-Te,0] [-rso] 
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where the constant Mg = P(E e (0)), i.e. the constant Mq only depends on initial data 
(|6.1cl) . By Lemma 13.21 and ()4.5[) , we see that Mq is bounded by a constant Mq which is 
independent of e, so that sup t6 r j. e i E e (t) < 2Mq. 
We therefore see that if we set 

T = mm (-*— °' + C+ M ° ^ (7 4) 

V 4Ci M ' 2Pi (2M ) ' 2P 2 (2M ) / ' 1 ' ' 



equation (17.ip implies that on [— T, 0], 
|r? £ (x,t) -rf(y,i)| > -\x-y\ for (a 
while equation (|7.2|) shows that on [— T, 0] 



^(^,0-^(2/,*)! > |k-2/l for (x,y) e 9f(B)x 01(B), (I,k)£{(-1,0),(0,-1)}, (7.5) 



|r ? e (x ) t)-7 ? e (y ! t)|>e+(C_ + C + )|- for all (x, y) G 61(B) x . (7.6) 

We then have from (|7.5p and (|7.6p that the domain rj £ (t,Q e ) does not self-intersect for 
each t € [— T, 0] and from (|7.4p we also have the estimate 

sup E e (t)<2M . (7.7) 
te[-T,o] 

Since T > is independent of e by (|7.4p . the estimates we have just obtained will permit 
the use of weak convergence to find the initial domain CIq at t = —T and the initial velocity 
field uq at t = —T, from which the free surface Euler equations, when run forward in time 
from t = 0, will produce the self-intersecting splash domain £l s and velocity field u s at the 
final time T > 0. 

8. ASYMPTOTICS AS € -> ON THE TIME-INTERVAL [— T, 0] 

8.1. Construction of the initial domain the asymptotic domain at t = —T. 

Lemma 16. 1 1 together with the estimate (|7.7p shows that 

\\Cif(-T,6f)\\l 5 , B <^M . 

l=-l 

Weak compactness and Rellich's theorem provide the existence of a subsequence (which by 
abuse of notation we continue to denote by r/ e ) such that 

rf(-T, •) o Q\ -± Qi , as e -> , in H 4 - 5 (B,) , (8.1a) 

rf(-T, .) o 0f -> 6j , as e -> , in H 35 (B,) , (8.1b) 



where B<; = B C\ B(0, q) and <j is given in Definition | 

We now define fio as the union of the sets 6;(£> f ) (— 1 < I < L). Due to (|7.5p . (|7.6p and 
(|8.1bp . we have that 

(i,I,)€B t xB t ,((I ) ^{(-l I 0) ) (0-l)}) |0,(a;)-e fc (j/)| > l|e,(a;)-0 fc (y)|, (8.2) 
and from (72]) on [-T,0], 

V(s,y)€jB f XjB s , |6_(x)-e + (y)| > (C_ + C+)|, (8.3) 
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where 0_ = 0_i and 0_|_ = ©o- These inequalities show that the boundary of f^o does not 
self-intersect and that £Iq is locally on one side of its boundary. Furthermore, setting k = I 
in (18. 2p . we see that each smooth map 0/ is injective, and thus each Qi(B f ) is a domain, 
which implies that Oo is an open set of M 3 . 

Lemma 8.1. is a connected, H 4 5 -class domain, which is locally on one side of its bound- 
ary. 

Proof. Step 1. We begin by proving that Qq is connected. To this end, fix X and Y in 0,$ 
so that X E @i(B q ) and Y E 0j(£ ? ) (-1 < l,j < L). We let (x,y) E B q x B q be such that 
X = &i(x) and Y = &j(y), and we define 

X* = ri e (6t(x),-T)eTi e (n',-T) 
Y e = T, e (ef( y ),-T)€T, e (n e ,-T). 

For (3 > 0, we set 

np = {z e n e | dist(z, on e ) > p} . 

Then for (3 > small enough, we have that Q% is connected and X e and Y e are in r] e (£lp, —T). 

From (f8lbl) we infer that each rf{9l, — T) uniformly converges to 0/ in B q ; thus, for e > 
small enough, we find that 

1 e (^,-T)cJ] , (8.4a) 
I e e6,(B ? ), r £ e^). (8.4b) 

Now, as dp is a connected set, so is n e (£lp, — T). Since X e and Y e are in this connected 
set, we let Cx e ,Y e denote a continuous path included in ^({Ip, —T), and having X e and Y e 
as its end-points. From (|8.4ap . Cx*,Y e c ^o- 

Next since both X and X e belong to the connected set 0/(£> f ), let Cx,x* denote a 
continuous path included in ®i{B<;) C fio an d having X and X e as end-points. Similarly, 
we let Cy*,y denote a continuous path included in 0j(£> f ) E i7o an d having Y e and Y as its 
end-points. We then see that the union of these three paths joins X to Y and is contained 
in Qq, which shows that fio is connected. 

Step 2. The fact that Qq is an ff 4 ' 5 -class domain follows immediately from the convergence 
given in (18. laH . 

Step 3. We conclude by showing that fio is locally on one side of its boundary, and that 
with 5° = J3°nS(0,?), 

dn = u£_ie,(s?) , (8.5) 

which will indeed complete the proof that Qq is a standard ff 4 ' 5 -class domain. 

To this end we first notice from (|8.1bj) and the fact that if is volume preserving, that for 
each I, 

detDQi = detD0j > Q > . (8.6) 

Also, from (|8.2p used when k = I, we notice that each 0^ is an injective map, which with 
(|8.6p provides 

d[Q l (B<)] = Q l (dB<). 

Therefore, 

dn e ufL.&idBj . (8.7) 
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Now, let us fix x G dB q n {x^ > 0}. We then have (since the only modified charts are 
modified close to the origin) that for any — 1 < I < K, 

ef(x) = e l (x). 

We also notice that there exists — l<k<L and y G Z3 ? such that 9i(x) = 9 k (y) (since 

o = u£ = _A(sj)- 

We also have that 9%{y) = 9 k (y), for otherwise k would be equal to either —1 or 0, in 
which case 9 k (y) would be in a very small neighborhood of xo, which, in turn, would imply 
that I must be equal to k (since the charts or 9+ do not intersect the other charts in a 
small neighborhood of xo), but then we would not be able to have x at a distance <j from 
the origin. 

We then have r] e (—T,9l(x)) = r] e (—T,9^(y)) which with (I8.1bp implies that 

@i(x) = e k (y)e@ k (B q )cn - (8.8) 

We can prove the same inclusion in a similar way if x G dB f and K + 1 < I < L. With 
(jO) . this yields 

^oCU^e,^ ). (8.9) 

Now, for X G @i(B®), we have X = @i(x), with x G B®. Now for any y G B q such 
that there exists —1 < k < L satisfying X = 0fc(y), we see from (|8.3p that (l,k) ^ 
{(-1,0), (0,-1)}. Therefore, from (E2J) we have 61, (a:) = 9 k {y). Then, as 90 = U^^fc^), 
we see that ?/ G 5°. Therefore, we have proved that X does not belong to U^ = _ 1 0/ c (^ ? ) = 
Slo- Thus X G <90o, which establishes (|8.5p . 

Together with (|8.6p and f)8.3[) , this establishes that f^o is a smooth domain locally on one 
side of its boundary, and concludes the proof. □ 



8.2. Asymptotic velocity at — T in the limit e — > 0. From our equivalence Lemma 16. II 
and (H2D, 

EKv%-T,9ml,B<^, 
i=-i 

which shows the existence of a subsequence (which we continue to denote by the index e) 
such that 

v e {-, -T)o9f^Vi as e , in H A {B, ) , (8.10a) 
v e {-, -T) o 9\ Vi as e , in tf 3 (£ f ) . (8.10b) 

We now define no on f2o as follows: 

V/ G {—1,0, 1, 2, ...,L} , Uo (e,) = V; onB ? . (8.11) 

In order to justify the definition in (|8.1ip . we have to check that if 0/(x) = Qj(y), for 
x and y in £> ; , then V,(x) = V/(y). We first notice that if 0/(x) = &j(y), then by (J8.3P 
we have (l,k) £ {(-1, 0), (0, -1)}. From (|8.2p . we then infer that 9i(x) = 9j(y) and thus 
|0f(x) - 9 e Ay)\ < c e , with lim e ^o c t = 0. 
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This then, in turn, shows that 

\v e (9f(x), -T) - v e (9*(y), -T)\ < \Dv%, -T)\ Loa ^ { __ T)) \ei{x) - 

< c € \Dv € (-, -r)| iO0 (ne(_T)) , 

which, thanks to (17. 7|) . implies that 

\v<(6f{x), -T) - v^e^y), -T)\ < c e M . 
By using (I8.10b|) . this then implies at the limit e — > 0: 

\Vi(x)-Vj{y)\<0, 

which concludes the proof. Also from (|8.1ip . we have that uq E H 4 (Qq), with ||ito||4,n — 
M . 

8.3. Asymptotic domain and velocity on (— T, 0] in the limit e — > 0. Having estab- 
lished the asymptotic limit as e — > when t = —T, we next consider the time interval 
(— T, 0). We employ the identical argument for taking the limit as e — > for the case that 
— T < t < as for the case that t = — T, leading to an asymptotic domain Q(t) of class 
i^ 4,5 and an Eulerian velocity field u(-,t) G H 4 (Q(t)) with \\u(t, OIUjQft) < Mo. 

At time t = 0, there is a slight difference in the asymptotic limit e — > 0, in the sense 
that the limit domain is the splash domain Q s , which is a self-intersecting generalized H - 
domain, with the corresponding limit velocity field is u s G -£f 4 ' 5 (O s ). This limit simply 
comes from the fact that \\6f — 0i\\i.5 b+ ~^ and ||it| o Q\ — u s o H4.5 s+ -> as e -> 0. 

8.4. Asymptotic Euler equations. It remains for us to prove that 

u f (t,x) = u(t-T,x), 0<t<T 
is indeed a solution of the free-surface Euler equations on the moving domain 

Q f (t) = fi(t - T) , 

which evolves the initial velocity uq and initial domain onto the final data at time t = T 
given by u s and Sl s . This will, in turn, establish the fact that after a finite time T, the 
free-surface of the 3-D Euler equations develops a splash singularity. 

We again consider the asymptotic limit as e — > 0. For each e > fixed, we solve the Euler 
equations forward-in-time using as initial data, £l e (—T) for the domain, and T) for 

the initial velocity. 

To this end, we first define the forward in time quantities for < t < T by 

n e f (t) = n e (t - t) , 

and 





= U e (;t 


-T) 






in n e f (t) , 


Vf(;t) 


= V e (;t 


-T)o 


V e (; 




in Q. £ f (0), 


v}(;t) 


= v e (-,t 


-T)o 




-T)- 1 


in O^(0), 


V%;t) 


= P e (;t 


-T) 






in Q e f (t) , 


q)(;t) 


= q e (;t 


-T)o 






in O^(0). 
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It follows that 

divn)=0 inf^(i), 

v j = u e f o rff = dtrfj in O)(0) , 

V}(;0)=e in O^(0). 

Prom the definitions of v £ , rj e , and u e in (|4.9p - (|4.1ip and by uniqueness of solutions to 
(|4.8p . we see that (v.pp'j) is a solution of (jl.ip on [0, T] with initial domain 0^.(0) and initial 
velocity «^(0), with the domain and velocity at time t = T equal to O e and u e s , respectively. 

In order to analyze the limiting behavior of these solutions as e — > 0, we write the Euler 
equations in Lagrangian form on the fixed domain B q by pulling back the equations from 
the reference domain 0^(0) using the following local coordinate charts: 

0f = 77 e (-r,0f) for I = -1,0, 1,2, 
Denoting the local inverse-deformation tensor by 

bt = [D(ri}oet)]- 1 , 
for — 1 < I < K, solutions of the Euler equations satisfy 

Jo 

d t v}oet + [bt} T D(q}oef) = 

q) o §f = 
fa},i$)o0f = (e,4(O))o0f 

together with 

^(T,n>(o)) = n 6 . 

For I = K + 1,...,L the same equations are satisfied with the exception of the boundary 
condition (J8J3H). 

Our a priori estimate Theorem I A. II shows that for each I = — 1, 0, 1, 2, L 

SUp (U f (t) O 9f\\l B + o + || g e (f) ^||2 A < ^ 

te[o,T] v ' 7 

where Mq is a constant that depends on the i/ 4,5 -norms of 9f and the ff 4 -norm of u e j(0). 

Thanks to Lemma 18,11 and the convergence in ()8.10|) . we see that Mq is bounded by a 
constant which is independent of e. As such, we have the following convergence in two weak 
topologies and one strong topology: 

vjoOf-^VfoQt, in L 2 (0,T;ii~ 4 (S f )), (8.15a) 

rf f o 9\ -)•)?/ o 9; , in L 2 (0, T; H 3 (B q )) , (8.15b) 

Qf ° Of q.f o 6 ; , in L 2 (0, T; F 4 ' 5 (^)) , (8.15c) 



in B ? x (0, T] , 


(8.14a) 


in B c x (0, T) , 


(8.14b) 


in B, x (0, T) , 


(8.14c) 


on B x (0,T), 


(8.14d) 


on B, x {t = 0} , 


(8.14e) 




(|8.l4n 
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o 



which together with the convergence in (f8lbl) shows that for I = -1, 0, 1, 2, K, the limit 
as e — > of the sequence of solutions to (I8.14p is indeed a solution of 

Q l = @ t + [ vjoSi in B s x (0, T] , (8.16a) 

o 6, + [bif D(q f o Oi) = in B f x (0, T) , (8.16b) 

div n/0 e, v f o 6/ = in B f x (0, T) , (8.16c) 

q f o6| = on S x (0, T) , (8.16d) 

fa/, «,) o 9; = (e, «q) o G ; on B, x {t = 0} , (8.16e) 

j7/(r,n ) = n., (8.i6f) 

where b/ = [D(r]f oQ;)]" 1 , and where uj, and ?yj are the forward in time velocity, pressure 
and displacement fields. 

A similar system holds for the interior charts 0/, with K+ 1 < I < L, with the exception 
of the boundary condition (I8.16dj) . Therefore, since the charts Qi define CIq, we have 
established that 

rj f = e+ f v f mn x(0,T], (8.17a) 
Jo 

dtvf + A T f Dq f = infi x(0,r), (8.17b) 

div Vf v f = in fi x (0, T) , (8.17c) 

q f = on <9S1 x (0, T) , (8.17d) 

(rjf,v f ) = (e,u ) on O x {£ = ()}, (8.17e) 
17/(^00) = ^, (8.17f) 

where the matrix ylj = [Z??7f] -1 . By a return to Eulerian variables this means that (uf,pj) 
is solution of (jl.ip with initial domain and velocity Qq and uq, respectively, and final domain 
and velocity at time t = T equal to the splash domain f2 s and u s . 

9. The splat domain fi s and its approximation by standard domains fi e 

9.1. The splat domain. Whereas our splash domain has a boundary which self- intersects 
a point xo, an obvious generalization allows to define the so-called splat domain fi s , with 
boundary <9fi s which self-intersects on an open subset To of dft s . 

9.1.1. The definition of the splat domain. 

(1) We suppose that To C T s := dfl s is the unique boundary self-intersection surface, 
i.e., fi s is locally on each side of To for each xq £Tq. For all other boundary points, 
the domain is locally on one side of its boundary. We assume the existence of a 
smooth level set function <j) G -ff~ 4 ' 5 (M 3 ) such that Tq C {4> = 0} 

(2) We let Uq denote an open neighborhood of xq in M 3 , and then choose an additional 
L open sets {U{\f =1 such that the collection {t/z}/£ is an open cover of T s , and 
{Ul}f= is an open cover of Q s and such that there exists a sufficiently small open 
subset wCf/o containing To with the property that 

uinUi = for all 1 = 1,...,L. 



24 DANIEL COUTAND AND STEVE SHKOLLER 

We set 

£7+ = u n n s n {4> > 0} and Uq = u n n s n {^< 0} . 

Additionally, we assume that Uq n fi s Pi {0 = 0} = To, which implies in particular 
that Uit and £7,7 are connected. 




Figure 5. Splat domain fl s , and the collection of open set 
{Uq, Ui, Ui, Uk} covering T s . 

(3) We furthermore assume that our level set function is such that ||-D0|| > Cq > on 
Uq. 

(4) For each / G {1, ...,K}, there exists an 7f 4,5 -class diffeomorphism 0/ satisfying 

0, : := B(0, 1) -> E/i 

f/inn s = fl,(B + ) and ^nr 8 = e,(s ), 

where 

S + = {(xi,x 2 ,x 3 ) G 5 : x 3 > 0} , 
B° = {( Xl ,x 2 ,x 3 ) eB : x 3 = 0} . 

(5) For L > K, let denote a family of open sets contained in Q s such that 
{Ui}^ =0 is an open cover of f2 s , and for I G {i7 + 1, #z : B — > U\ is an i7 45 
diffeormorphism. 

(6) To the open set Uq we associate two i7 4,5 -class diffeomorphisms 9 + and of S 
onto Uq with the following properties: 

6 + (B + ) = U+, 9^B+) = Uq, 

o + (B°) = u+nr s , e^B ) = u z nr s , 
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such that 

T o = e + (B°)n0^(B o ), 

and 

+ = 0_ on uJo c -Bo • 
where uq is a smooth connected domain of Bq in M 2 . 
We further assume that 

6±(B+ n 5(0, 1/2)) n 0jCB+) = for Z = 1, 

and 

6±{B+ n B(0, 1/2)) n 9jB) = for Z = K + 1,...,L. 

Definition 9.1 (Splat domain fi s ). We say that fi s is a splat domain, if it is defined by 
a collection of open covers {Ui}fL and associated maps {9±, 9±, 02, 0l} satisfying the 
properties (l)-(6) above. Because each of the maps is an H 45 diffeomorphism, we say that 
the splat domain fi s defines a self-intersecting generalized H 4 5 -domain. 

9.2. A sequence of standard domains approximating the splat domain. We ap- 
proximate the two distinguished charts 0_ and 0+ by charts 01 and 6+ in such a way as to 
ensure that 

e e _(B°)ne e + (B°) = v e > 0, 

and which satisfy 

61 -»• 0_ and 9% -»• 0+ as e -> . 

We let ip G ^(w) denote a smooth bump-function satisfying < ip < 1 and i/) = 1 on To- 
For e > taken small enough, we define the following diffeomorphisms 

et( x ) = 0_(x) - e ij>(0-(x)) D<j>{6-(x)) , 

9%{x) = 0+(x) + e il>(d+{x)) D<p(9 + (x)) , 

By choosing ip G T>(lj), we ensure that the modification of the domain is localized to a small 
neighborhood of To and away from the boundary of U$ and the image of the other maps 9i . 
Then, for e > sufficiently small, thanks to item (3) in the definition of the splat domain, 

(b(8L(x)) < 0(0-0*0) " f ^(MaO) WU*))| 2 < 0, 

i>(e%(x)) > H9 + (x)) + e -^(6+(x)) \D<i>{9 + {x))\ 2 > , 

which shows that 

9L{B+) n 9%{B+) = 0. 

Since the maps 0|_ are a modification of the maps 9± in a very small neighborhood of G B, 
we have that for e > sufficiently small, 

0±(5 + n B(0, 1/2)) n 9i(B + ) = for Z = 1, if , 

and 

9 e ±(B + n 5(0, 1/2)) n 9i{B) = for Z = if + 1, L . 
For Z G {1, we set 0f = 0;. Then 01 : B + -> £/ , 0+ : 5+ -> Z7 , and 0f : 5+ -> tf,, 

Z G {1, ...,-ftT}, 0[ : B — > Ui, I G {K + 1, is a collection of H 45 coordinate charts as 

given in Section 12.51 and so we have the following 
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Lemma 9.1 (The approximate domains fi e ). For each e > sufficiently small, the set fi e , 
defined by the local charts 01 : B + -> Uq, 9% : B + -> C/ , arid 0f : B+ ^ U h I £ {1, if}, 
9f : B Ui, I £ {K + 1, L} (given in Definition \9.1\) is a domain of class H 5 , which is 
locally on one side of its H 4 boundary. 

Just as for the splash domain, we have approximated the self-intersecting spZai domain 
fi s with a sequence of ff 4 ' 5 -class standard domains fi e locally on one side of its boundary 
for each e > 0. Also, just as for the splash domain, our approximate domains fi e differ from 
our splat domain f2 s only on the two patches 9 t _{B + ) and 9+(B + ). In particular, as 9± 
differ from 9± on a set properly contained in uj C Uq, we continue to use the same covering 
{Ui}{l q for fl e as for $l s . 

10. Construction of the splat velocity field u s at the time of the splat 

singularity 

We can now define the splat velocity u s associated with the generalized H 4,5 -class splat 
domain f2 s , as well as a sequence of approximations u| set on our i? 4 ' 5 -class approximations 
fi e of the splat domain Q s . 

10.1. The splat velocity u s . 

Definition 10.1 (Splat velocity u s ). A velocity field u s on an H 4,5 -class splat domain Cl s 
is called a splat velocity if it satisfies the following properties: 

(1) Cu s o 9± <E H 4 - 5 (B + ), (u s o 6i G H 4 - 5 (B+) for each 1 < I < K and u s G F 4 - 5 (w) /or 
eac/i aJ C fi s ; 

f^J u s • Dcj)o 0_|e_(B+) > C- ^d — u s • L>0o 0+|e + (B+) > C + with C- + C + > 0, so that 
under the motion of the fluid, the sets Uq and Uq are moving relatively towards 
each other. 

We can then define the approximate splat velocity fields u| : Og — > M 3 in the same way as 
we did for the case of the splash velocity. The results of Sections [7] and [8] can then proceed 
in the same fashion as for the splash case, leading to Theorem 15.21 

We note only that the inequality ()7.2p must replaced with 

- ^(x, t)) + <f,(r, e (y, t)) > (C+ + C_) \t\ - t 2 P(sup E") , (10.1) 

[0,t] 

for x, y as in (|7.2p . The estimate (jlO.ip together with 

IWOM)) - Hv e (y,t))\ < \D<p\ \ V e (x,t) - V e (y,t)\ , 
and item (3) of the definition of our splat domain fi s then provides 

\rf{x,t)-tf{v,t)\ > {C+ + C - ] |t| -t 2 Pi(sup^). (10.2) 

G o [0,t] 

In turn, (|10.2p allows to establish e-independent estimates and arrive to the analogous 
conclusions as those obtained in Sections [7] and [8l 
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Appendix A. A priori estimates for the free-surface Euler equations 

In this appendix, we establish a priori estimates for the free-surface Euler equations with 
reference (or initial) domain Q which is a standard i/ 4,5 -class domain, open, bounded, and 
locally on one side of its boundary. 

A.l. Properties of the cofactor matrix a, and a polynomial-type inequality. 

A. 1.1. Geometry of the moving surface T(t). With respect to local coordinate charts, the 
vectors r), a for a = 1,2 span the tangent space to the moving surface F(t) = rj(F) in IR 3 . 
The (induced) surface metric g on T(t) has components g a p = r], a -T),p. We let go denote 
the surface metric of the initial surface T. The components of the inverse metric are 
denoted by [g] al3 ■ We use y/g to denote \/det g; we note that yfg = \r\,\ xrj,2 |, so that 
n(v) = x V,2]/x/g- 

A. 1.2. Differentiating the inverse matrix A. Using that Dr\A = Id, we have the following 
identities 

dA k = -Aid v r , s A k , (A.l) 
DA* = -AfDri r , s A k r , (A.2) 
d t A k = -Afv r , s A k . (A.3) 

A. 1.3. Relating the cofactor matrix and the unit normal n(t). With N denoting the outward 
unit normal to T, we have the identity 

ni ( V ) = a k N k /\a T N\ . 

so that 

AfN k = J^yfgmifi) onT. (A.4) 

A. 1.4. A polynomial-type inequality. For a constant Mq > 0, suppose that f(t) > 0, t t— > 
f(t) is continuous, and 

f(t)<M + tP(f(t)), (A.5) 

where P denotes a polynomial function. Then for t taken sufficiently small, we have the 
bound 

fit) < 2M . 

We use this type of inequality (see [13]) in place of nonlinear Gronwall-type of inequalities. 

A.2. Trace and elliptic estimates for vector fields. The normal trace theorem states 
that the existence of the normal trace w • N\r of a velocity field w G L 2 (J7) relies on the 
regularity of divw G L 2 (Q) (see, for example, [24] )• If divw G L 2 (Q), then w • exists in 
i7" 5 (r). We will use the following variant: 

\Bw ■ Af_ a5ir < c[\\dw\\l n + lldivHIIU] ( A - 6 ) 

for some constant C independent of w. 

The construction of our higher-order energy function is based on the following Hodge-type 
elliptic estimate: 
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Proposition A.l. For an H r domain 0, with T = dQ, r > 3, if F G L 2 (Q;R 3 ) with 
curlF e H^^R 3 ), divF G H 8 ' 1 ^), and dF ■ N\ T € H S ^(T) for 1 < s < r, then there 
exists a constant C > depending only on f2 such that 

\\F\\ s ,a < C (\\F\\ ,n + \\ cur\F\\ s _ hn + || divF|| s _ lin + \dF ■ AVf.r) . ( A - 7 ) 

where N denotes the outward unit-normal to T. 

This well-known inequality follows from the identity —AF = curlcurli 7 — DdivF. 

A. 3. The higher-order energy function E{i). 

Definition A.l. We set on [0, T] 

E(t) = 1 + h{t)\\l^ n + \\v(t)\\l,n + II ^<t)\\i. 5 ,n + lh(*)lli.5,n • (A-8) 

The function E(t) is the higher-order energy function which we will prove remains bounded 
on [0,T]. 

Definition A. 2. We set the constant Mq to be a particular polynomial function P of E(0) 
so that M = P(E(0)). 

A. 3.1. Conventions about constants. We take T > sufficiently small so that, using the 
fundamental theorem of calculus, for constants c\,C2 and t £ [0,T], 

-{qo9i), 3 (i)>A/2 Vl = l,...,K, 

ci det g(0) < det g(t) < c 2 det g(0) on T , 

||r?(£)|| 4 < l|e|| 4 + l, ||g(t)|| 4 < |k(0)|| 4 + l, 

\\v(t)\\ 3 .5 < ||«o||3.5 + 1 , \\vt(t)\\ 3 < \\vt(0)\\ 3 + 1 • 

The right-hand sides appearing in the last three inequalities shall be denoted by a generic 
constant C in the estimates that we will perform. The norms are over O. 

A. 4. Curl and divergence estimates for w, v, and Vf 

Proposition A. 2. For all t e (0,T), 

l|curlr/(t)|||. 5in + ||curl^(t)||l. 5|n < M + TP( sup E(t)) . (A.9) 

te[o,T] 

Proof. By taking the curl of (|1.3b|) . we have that 

cur I,, Vt = . 

It follows that dticwclqv) = B(A,Dv), where the feth-component of B is given by 
[B(A, Dv)] k = EkjiAt'jV^s = e kij v\ s A s p v p d A\ ; 

hence, 

cm\v{t) = curlu + / B(A(t') , Dv(t'))dt' . (A.10) 
Jo 

Step 1. Estimate for curl?]. Computing the gradient of (jA.lOp yields 

cml n Dv(t) = DcutIuq - e. ji DA s j v i , s + DB(A(t'),Dv(t'))dt' . (A. 11) 
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(In components, [curl^ d Xl v]i = £ijkV k ,i r Aj.) Applying the fundamental theorem of calculus 
once again, shows that 

curl v Di](t) = tD curl uo + e.jt [ [A t jDrf, s —DA s ,v i , s ]dt' 

Jo 

+ f f DB{A(t"),Dv{t"))dt"dt' , 
Jo Jo 

and finally that 

D curl r](t) = tD curl uq — e.ji f A t s At')dt' Drf, a (A.12) 

Jo 

ft ft ft' 

+ e.ji / [Atprf^-DAy^dt + / / DB(A(t"),Dv(t"))dt"dt' . 
Jo Jo Jo 

Using the fact that dfA? = —Afv l , p A P and DA? = —AfDr] l , p A p , we see that 

DB(A,Dv) = -e kji [Dv\ a Atf, p A p j +v i , 8 AfDj, p A* 

+ v i , s v l , p D(AfA p j )]. 

The precise structure of the right-hand side is not very important; rather, the derivative 
count is the focus, and as such we write 

DB(A, Dv) ~ D 2 v Dv A A + D 2 r] DvDvAA. 

Integrating by parts in time in the last term of the right-hand side of (IA.12P , we see that 

i-t rt' ft ft' ft ft' 

I \ DB(A,Dv)dt"dt' ~ - / / D 2 r](DvAA) t dt"dt' + / D 2 r] Dv Dv AAdfdt' 
Jo Jo Jo Jo Jo Jo 

+ / D 2 rj Dv A Adt' . 
Jo 

Thus, we can write 

D curl rj(t) ~ tD curl u + D 2 r] f Dv A Adt' + f D 2 rjDvAAdt' 

Jo Jo 

ft rt' 

D 2 r]DvDv AAdfdt' + / / D 2 r] (Dv AA) t dt"dt' . 

Jo Jo 



JO 

Our goal is to estimate H^D curl r7||| 5 ^- Thanks to the Sobolev embedding theorem, we 
have that 

\\D cuvl v \\l^ u < M +TP{ sup E(t)), 

te[o,T] 

and hence with curl^ vt = 0, that 

l|curl7/|||. 6)n <Mo + rP( sup E(t)). 

te[o,T] 
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Step 2. Estimate for curL v. Integrating-by-parts with respect to dt in the time in- 
tegral in equation (lA.lljl , we see that the highest order term in curb, Dv is given by 
£D 2 v Dv t AAdt'. As H 2 - 5 (n) is a multiplicative algebra, it follows that on [0, T], 

||curl^(t)||§. 5 . n <M + TP( sup E(t)) . 

te[o,T] 

□ 

Proposition A. 3. For all t e (0,T), 

\\dw V (t)\\l 5>n + \\dwv(t)\\in<M + TP( sup E{t)). (A.13) 

te[o,T] 

Proof. Since A\v % ,j = 0, we see that 

AlDv i , j =-DA>v i , :j . (A.14) 
Step 1. Estimate for divry. It follows that 

;,H/V.,/ ih^Dn,, DAlr'., . 
Using the fact that r](x, 0) = x, 

[AiDn^jKt) = £ {ihA\Dn'.j -da{ v\ 3 ) dt' , 

and hence 

Ddwri(t)= [ d t A{Dr)\jdt' - I DA\v\idt! - [ d t A{dt'Drf, 
Jo Jo Jo 

Again, the Sobolev embedding theorem provides us with the estimate 

\\diyr,(t)\\i 5jQ <TP( sup E(t)) . 

te[o,T] 

Step 2. Estimate for divv. Prom A\v l ,j = 0, we see that 

divv(t) = -[ d t A\di! v\j . (A.15) 
Jo 



j ■ 



Hence 

\\diw(t)\\l n <TP( sup E{t)). 

te[o,T] 



□ 



A. 5. Pressure estimates. Letting A\-J^- act on (jl.3b|) . for t G [0,T], the Lagrangian 
pressure function q(x,t) satisfies the elliptic equation 



Ah, 



,j = v\j A{ v r , s Af in n , (A.16a) 
q = on T. (A.16b) 
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By invoking the Sobolev embedding theorem, elliptic estimates in our local charts show 
(with the norms set on Q) 

\\q\U<C(\\A\\ 2 ,\\vh)\\A\\ 3 , 
\\qh<C(\\A\\ 2 ,\\v\\ 3 )\\A\\ 4 , 

where the constant has polynomial dependence on ||j4||2 and \\v Linear interpolation 
then yields 

IMI4.5 < C(PI|2, 1Mb) IMU.5 • 
By time-differentiating (|A.16p , and using our conventions of Section IA.3.11 concerning the 
generic constant C, we have the elliptic estimate on [0, T] 

||?(t)||4.5 + ||ft(t)||4<C||T;(t)||4.5. (A.17) 

Remark 2. When the elliptic problem ( L4. 16]) is set on the approximate splash domain £l e , 
the elliptic constant a priori depends on e > 0, via the charts 9±; however, thanks to Lemma 
\3.2\ the elliptic constant is independent of e since the charts 9± are bounded in H . 



A. 6. Technical lemma. Our energy estimates require the use of the following 

LEMMA A.l. Let Hz(£l)' denote the dual space of H^iVL). There exists a positive constant 
C such that 

\\dF\\ Hh{a)i <C\\F\\ hn VFG^(O). 
Proof. Integrating by parts with respect to the tangential derivative yields for all G £ 



fdFGdx = y"f Q[(F oOfi^oOi 1 Gdx 
Jn l=l Juinn 

L 

= J2 C°0i(Fo0i),a Go Oldest DOidx 
i=i Jb+ 

L 

= -J2 C°0lFo0i(Goei) >ol detDeida 
i=i Jb+ 

-J2 f FoeiGoei^oQid^De^dx 
i=i Jb+ 

<C||F|| 0l n||G||i ) n J 
which shows that there exists C > such that 

VF£L 2 (0), \\dF\\ H i {n y <C\\F\\ 0>Q . (A.18) 
Interpolating with the obvious inequality 

yF€H\n), \\BF\\ L 2 {n) <C\\F\\ 1>Q 
proves the lemma. □ 
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A. 7. Energy estimates for the normal trace of i] and v. By denoting rji = rj o 9\ we 

see that 

r]i(t) : B + -> n(t) for Z = 1, 

We set u; = uorji, qi = po rji and = [Dry/] -1 , J; = detDr/i, and a; = JiA\. It follows that 
for I = 1,...,K, 

Vi(t) = 0i+ [ vi mB + x[0,T], (A.19a) 
J o 

+ = in B + x (0, T] , (A. 19b) 

div^ uj = inB + x [0, T] , (A. 19c) 

q t = onfi°x[0,T], (A.19d) 

(Vi,vi) = (e l ,u o9 l ) inB + x{t = 0}. (A.19e) 
Proposition A. 4. Fort e [0,T], 

|ar ? (t)-iV|l r + |^(t)-iV|| 5r <Mo + rP( sup E(t)). (A.20) 

*e[o,T] 

Proof. We compute the following L 2 (B + ) inner-product: 

= (Cd^dm + A, D qi ] , (B\) L2{B+) . 
To simplify the notation, we fix I € {1, ■ ■■,!(} and drop the subscript. We have that 

= ~HCaV*)lloi* + / ^A k q, k B A v l dx + [ ( 2 A k B 4 q, k B 4 v l dx+TZ, (A.21) 

V v „ ' > v ' 

11 x 2 is 

where 7Z denotes integrals over B + consisting of lower-order terms (or remainders) which 
can easily be shown, via the Cauchy-Schwarz inequality, to satisfy 

cT 

\K{t)\dt < M + TP{ sup E(t)) . 

te[o,T] 

Using the identity (|A.1|) . we see that 

^A^BS^sAfq^B^dx + n 

( 2 A k B* V r q, k B A v i A s i N s dx h + [ CA k r BS r A\q, k 0V„ dx +1Z , 

B° JB+ 



where dx^ = dx\dx2 denotes the surface measure on B°. As q = on B°, q a = and 
q,2 = on T, so that 

l2a = I (~q,3 )C 2 B 4 r, r A k r N k B A v'Af dx h . 
Jb° 

We define n\ to be the outward unit normal to the moving surface rji (t, B°), so that from 

dA3D, 

AitN k = Jf 1 ^gin l (r H ) on B° . 
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Dropping the subscript I again, it follows that 

Z2a(t) = / (-<7,3 )C 2 <9 4 ?? • ™ d A v -n\ det g\J~ 2 dx h 
Jb° 

= f (-9.3 K 2 |9 4 7? • n| 2 | det g\J- 2 dx h - [ VdVdy | det g\ J~ 2 ]dx h . 

l at j B o J B o i 

* . ' > , ' 

K a K-b 

By the assumption of Section IA.3.11 

\dt\nirij \ det g\J- 2 ]\ Loc ^ r) <C, 



from which it follows that 

r-T 



I K b {t)dt <CTP{ sup E(t)). 
Jo telo.T] 



'o te[o,T] 
Using our assumed bounds for —q,s (t), det g(t), J on [0, T], we see that 

c\(B\(t) ■ ni (t)\ 2 B0 < [ x 2a (t)dt, 

Jo 

for a constant c which depends on A,g(0), and J(O) = det D9i. We set 

jy _ Phi X 0l,2 

By the fundamental theorem of calculus ni(t) = Ni + Jq dtni(t')dt' , and by our assump- 
tions in Section lA.3.11 supj 0T ] \9tni(t)\ Lao ^ < C; hence, 



c\(d\(t) ■ N^l Bi) < ! l 2a (t)dt + TP( sup E(t)), 
Jo te[o,T] 



and hence 

r-T 



c|C^l(Q-M|l B 0< / l 2 a(t)dt + TP{ SUp 

Jo tefo.n 



'0 t€[0,T] 



It remains to show that the integrals f Q X2b{t)dt and J Q X^{t)dt are both bounded by 
TP(su Pte[0iT] E(t)). Using (|XT9cl), 

X 26 (t) = - / C 2 ^ #Y u<„ d 4 ,4| cte + K 
Jb+ 

<c\\cdS(t)\\ hB+ \\<8 4 Am H i {B+Y +^ 

<c\\d 4 v (t)h-\\d 3 A(t)\\ H , {n)+ n 

< C sup E(t) + TZ , 
te[o,T] 

where we have used Lemma IA.ll for the second inequality. 
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Finally, 

X 3 (*) = -/ C 2 d A q B A v\ k A\ dx = [ d A qv i , k d A A k i dx + n 
Jb+ Jb+ 

<c||ca 3 ^)lli,B + IIC9 4 A(t)|| H i (B+)/ +^ 

< C sup E(t) + K , 

te[o,T] 

where we have used the pressure estimate (|A.17|) and Lemma lA.ll for the last inequality. 

Summing the estimates for I±, I2, I3 and integrating (|A.21|) from to T, we obtain the 
inequality, 

sup (\(dv(t) ■ ml B0 + \\Cd A v(t)\\l B+ ) <M + TP( sup E(t)). 

te[o,T] v ' 7 *e[o,T] 

According to Proposition IA.31 



sup ||divu(t)||§ < M + TP( sup E(t)) 

t£[0,T] te[0,T] 



from which it follows that 



sup \\(d 3 divv(t)\\l B+ <M + TP( sup E(t)) 
te[o,T] te[o,T] 



Hence, the normal trace estimate ()A.6P shows that 



sup ICSM^-Mliiflo <M + TP(sup 
te[o,T] v 2 ' 7 te[o,T] 



from which it follows that 



sup (\Cdri(t)-Ni\l^ + \Cdv(t)-Ni\l BBO )<Mo+TP(8vip E(t)) . 
*e[o,T] v ' 7 te[o,T] 



□ 



Combining Proposition IA. 41 with the curl estimates in Proposition I A.2I and the divergence 
estimates in Proposition IA.3I for rj(t) and v(t) and using ()A.7|) together with the fact that 
vt = —A T Dq provides us with the following 

Theorem A.l. Suppose that the initial pressure pq satisfies < on T and that E(0) < 
00. For T taken sufficiently small and for a polynomial function Pi, 

sup (h(t)||!.5,Q + lbWllln + l|curl^(t)||| 5in + || Wt (t)||| 5in ) <M + TP 2 ( sup E(t)) . 



te[o,T] te[o,T] 
Moreover §E < on r(t) /or t G [0,T]. 
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